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^ ■ Abstract 

' , We investigate whether the universe was homogeneously in the false vac- 

D , uum state at the critical temperature of a weakly first-order phase transition 

' such as the electroweak phase transition in terms of a series of numerical simu- 

lations of a phenomenological Langevin equation, whose noise term is derived 
^ ■ from the effective action but the dissipative term is set so that the fluctuation- 

dissipation relation is met. The correlation function of the noise terms given 
by a non-equilibrium field theory has a distinct feature if it originates from 
interactions with a boson or with a fermion. The spatial correlation function 
of noises from a massless boson damps with a power-law, while the fermionic 
noises always damp exponentially above the inverse-temperature scale. In 
the simulation with one-loop effective potential of the Higgs field, the latter 
turns out to be more effective to disturb the homogeneous field configuration. 
Since noises of the both types are present in the electroweak phase transition, 
our results suggest that conventional picture of a phase transition, namely, 
nucleation of critical bubbles in a homogeneous background does not apply 
or the one-loop approximation breaks down in the standard model. 



I. INTRODUCTION 



In modern cosmology our universe has presumably experienced a number of phase tran- 
sitions in the early stage of its evolution. Among them those at the grand unification scale, 
which may be related with inflation and/or formation of topological defects, are very spec- 
ulative in that we know neither the symmetry-breaking pattern nor the initial state before 
the transition which must be highly non-thermal due to the rapid cosmic expansion then. In 
contrast, the dynamics of the electroweak phase transition (EWPT) in the standard model 
is a much more solid subject of study because we have a perturbative expression of the 
effective potential with only one undetermined degree of freedom, namely, the Higgs mass, 
Mhi [0 01 ^'^d ^6 may well assume the thermal state before the transition since the 
relevant particle interaction rates are much larger than the cosmic expansion rate by that 
time. 

Nevertheless the dynamics of EWPT is not fully understood yet, mainly because, al- 
though one-loop effective potential of the Higgs field, 0, shows it is of first order, the po- 
tential barrier between the two minima at the critical temperature, Tc, is so shallow that it 
has been doubted if the conventional picture of nucleation of critical bubbles in the homo- 
geneous false-vacuum background really works. Whether the transition is of first order with 
super-cooling is a very important cosmological issue to judge if electroweak baryogenesis 
is possible. 

Much work has already been done on this topic. For example, Gleiser, Kolb, and Watkins 
^ considered the role of subcritical bubbles of the correlation volume as a noise effect in a 
weakly first-order phase transition, and Gleiser and Kolb concluded that for the Higgs 
mass larger than 57GeV, the universe is not in the false vacuum state uniformly with = 
but in a mixture of = and the true vacuum = 0+(Tc)|^. See also Gleiser and Ramos 
0. Gleiser and Borrill and Gleiser ^ have confirmed occurrence of such "phase mixing" 
by numerical simulation of a phase transition, solving a simple phenomenological Langevin 
equation of the form 

□0(a;,t) + VM.Tc) = i^x^t) , (1.1) 

at the critical temperature T = Tc. Here an overdot denotes time derivation, Vefr(0, Tc) is 
an effective potential, and C,{x,t) is a random Gaussian noise with the correlation function 

( ^{Xi, ti)^{X2, h) ) = V6{Xi - X2)6{h - t2) , (1.2) 



satisfying the fluctuation-dissipation relation V = 2riTc. Later Shiromizu et al. ]T0|| treated 



the size of a subcritical bubble as a statistical variable and discussed its typical size is 
smaller than the correlation length. They concluded that for any experimentally allowed 
value of the Higgs mass, or Mh > 60GeV []TT|, the phase mixing does occur already at the 



*At the critical temperature, the notion of the false and the true vacua are not well defined since 
the both minima are degenerate. Here we use the same terminology as in lower temperatures. 
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critical temperature. Furthermore, in the Monte-Carlo lattice simulations of the Euclidean 
four dimensional theory or the reduced three dimensional model of the finite-temperature 
electroweak theory, analytical cross-over behavior is observed for most values of the Higgs 
mass permitted from the experiment [|12[ |13 . 



On the other hand. Dine et al. calculated root-mean-square amplitude of the Higgs 
field on the correlation scale, or the inverse-mass length scale at = at the critical temper- 
ature and concluded that for the Higgs boson with Mh — 60GeV it is much smaller than the 
distance between the two minima and the fraction of the asymmetric phase of the Universe 
is negligible (e~^^) so that subcritical fluctuations do not affect the dynamics of EWPT. Bet- 
tencourt |T^ confirmed their conclusion by estimating the probability that the mean value 
of averaged over a correlation volume is larger than the distance to the maximum of the 
effective potential separating the two minima and finding that it is extremely small. Finally 
in response to Shiromizu et al. [0], Enqvist et al. treated both the amplitude and the 



spatial size of subcritical fluctuations as statistical variables and discussed that subcritical 
bubbles, if they exist at all, resemble the critical bubbles and that the usual description of 
a first-order phase transition applies. Their analysis, however, has a problem that it suffers 
from severe divergence and they had to introduce cut off ad hoc. 

Thus a number of independent analyses have drawn different conclusions about how 
EWPT proceeds. In the present paper we attempt to elucidate why such discrepancy has 
arisen with the help of numerical simulations of a phenomenological Langevin equation 
which is better motivated than (|1.1|) and ( |1.2|) from a non-equilibrium field theory. In fact 
the origin of the discrepancy is quite simple: it only reflects at which spatial scale one 
estimates the amplitude of fluctuations. We try to approach the problem step by step. 

The rest of the paper is organized as follows. We start with a re-analysis of the sim- 
ple Langevin equation ( |l.lj ) with white noise ( [1. 21 ) in Sec. First we consider a non- 
selfinteracting massive scalar model and show that numerical calculations of (|1.1| ) and ( |1.2|) 



on a lattice can reproduce the finite-temperature spatial correlation function of a massive 
scalar field correctly as long as we take the lattice spacing comfortably smaller than the 
correlation length or the Compton wave length. We then solve the same equation but with 
an effective potential in the standard model as was done by Borrill and Gleiser |p. We find 
not only the behavior of the correlation function but also the limit on Mh above which 
phase mixing occurs change drastically depending on the lattice spacing. In order to ob- 
tain a sensible bound on Mh, therefore, we should choose a reasonable value of the lattice 
spacing with the help of a fundamental theory. This issue is discussed in Sec. ^ and a 
new phenomenological Langevin equation is proposed. In Sec. |^ we report the results of 
numerical simulations of the dynamics of the field based on this equation. In Sec. |V] we give 
an intuitive explanation for the numerical results by using a simple Boltzmann equation. 
Finally Sec. ^ is devoted to summary and discussion. 

II. RE- ANALYSIS OF THE SIMPLE LANGEVIN EQUATION 



In this section we elucidate the origin of the discrepancy in the previous literatures. For 
this purpose, we first solve the simple Langevin equation (|1 . 1|) in the case only the mass 
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term is present in the potential, namely, 



[D + m^) (p{x, t) + r]^{x, t) = ^{x, t) 



(2.1) 



where C,{x,t) is a random Gaussian noise satisfying ( |1.2| ) with T> = 2riT. 

After discretizing the system on a lattice, we follow the time evolution from the initial 
condition, (f){x, 0) = and (f){x, 0) = 0, on each lattice point. The dimensionless variables 
X = Tx, t = Tt, $ = T^'^(f), fj = T~^r], ^ = T~^^, and fi = T'^m are introduced for 
numerical calculations, but we omit the tildes below. We arrange three different lattices 
for comparison. One has the total lattice points, = 32^, grid spacing, 6x = 1.0, time 
step, 6t = 0.1, and total run time, t = 500, another has = 40^, 6x = 0.8, 6t = 0.1, and 
t = 500, and the other has A^ = 64^, 6x = 0.5, 6t = 0.1, and t = 500. Using the second 
order staggered leapfrog method, the discretized master equation reads 



s=x,y,z 



t,n X,n+l/2 ' 



{6x) 



(2.2) 



where i represents spatial index and n temporal one, and /i is set to be 0.125. The correlation 
of the noise is given on the lattice by 



{6x[ 



(2.3) 



Since it is white both spatially and temporally, we have only to generate Gaussian white 
noise on each grid. 



e 



27] 



1,11 



5t{5x) 



(2.4) 



where ^ is a Gaussian random number with a vanishing mean and a unit dispersion. The 
periodic boundary condition is imposed. Under the above conditions we take the ensemble 
average over five different noises for each case. The correlation function C(r), = ( ^{x)^{y) ) 
with r = \x — y\, is numerically obtained by calculating all the combinations of the product 
^(x)^{y) satisfying r — 0.5 < |a:; — y| < r + 0.5. 

On the other hand, the analytic expression of the equal-time correlation function of a 
non-selfinteracting scalar field with mass m at temperature T is given by 



1 

m 



E 



(27r)3a;2 + fc2 + m2 



,ik-{x - y) 



47r^r 



Ki{mr) + 2Y, 



(2.5) 



where r = \x — y\^ 13 = 1/T, ujn = 2?™//?, and Kj is the modified Bessel function of the j-th 
order. This correlation function damps exponentially above the inverse-mass scale. 
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Both numerical and analytic results are depicted in Fig. |l] in the dimensionless unit [[ We 
find that the correlation functions obtained from numerical calculation damp in a manner 
independent of the lattice spacings and also that they coincide with the analytic formula 
( p.5|) , namely, exponentially damp with the correlation length of the inverse- mass. Thus 
for the case of the massive non-selfinteracting scalar field, not only the simple Langevin 



equation ( |1 . 1|) with the random noise (|1.2| ) but also its numerical solution on a lattice can 
reproduce its actual finite-temperature behavior. 

Next we consider an interacting scalar field borrowing the one-loop improved effective 
potential, Vew: of the Higgs field in the electroweak theory following Borrill and Gleiser [Q. 
The master equation is given by 



t) + r]<}){x, t) + TJ = ^{x, t) , 

where ^{x,t) is again a random Gaussian noise with no correlation. Vew is given in 
i as 



(2.6) 



V, 



EW\ 



(2.7) 



where 



^-24 



E 



1 



6 



a 

M, 



6( — 

o 



0.169, 



w 



a 



a 



0.00965 



(2.8) 
(2.9) 



for M, 



w 



80.6 GeV, Mz = 91.2 GeV, Mt 



174 GeV, and a = 246 GeV |11|. We also find 



To 



Ml 



AD 



Mjj = (2Ao + 12S)(T^ 



B 



-(6M^ + 3M| - 12M/ 



647rV4 

and the temperature-corrected Higgs self coupling is given by 



0.00456 



Ar — An 



1 



167r2 



Y.9B 



B 



Mf 



a 



\ ( Ml\ ^ 



cr 



Am 



(2.10) 

(2.11) 
(2.12) 

(2.13) 



where the sum is performed over bosons and fermions with their degrees of freedom gB(F) 
and IncB = 5.41, In Ci;' = 2.64. The Higgs field, appearing in the potential ( p.7|) in the 
actual electroweak theory is of course the amplitude of an SU{2) doublet complex scalar 
field. But in solving the Langevin equation (|2.6|), we neglect the gauge-nonsinglet nature 



^The correlation function derived analytically is reguralized by subtraction of the vacuum energy. 
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of the field for simplicity and constrain its dynamics along the real-neutral component to 
treat (j) as if it was a real singlet field as in Borrill and Gleiser |^. Thus this should not be 
regarded as the simulation of the actual Higgs field, although our results would be suggestive 
to it. 

Introducing dimensionless variables, x = {2DY^'^T2X, t = (2Dy^^T2t, $ = 
(2D)-V4T-V, V = {2D)-y^T2 \ I = {2D)-'/^T2% 6 ^ T/T^, and VEwi(l>) = 
{2D)'^/^TiU{^), Eqs. and Q reduce to 

0^x) + ^^^ + r]^x',i) =^{S:) , (2.14) 

and 

f/($) = 1(02_^)^2_|^^3^^^4_ ^2.15) 

Here dimensionless parameters are defined as a = (2D)~^^^{3E) = 0.065, A = {2D)^^^'^\t = 
1.72At. Hereafter we omit the tildes. The effective potential is depicted in Fig. For 
9 > 9i = {1 — a^/4A)~^/^ there is only one minimum at $ = 0. At 6' = appears the 
inflection point $ = a6i/2X. As the temperature further drops, another minimum, which 
is metastable, appears and at the critical temperature, 9c = (1 — 2a'^ /9\)~^^'^, the two 



minima, $ = 0, and $+ = || 1 + - 4A(1 - 1/9"^) /a^ , are degenerate. Below 9 = 9, 



'CI 



the symmetric state $ = becomes metastable in turn and at ^2 = 1 disappears the local 
maximum at <!>_ = || 1 — ^jx — 4A(1 — l/9'^)la'^ . This is a typical model which represents 
a first-order phase transition. 

We investigate whether the universe is in a homogeneous state of the false vacuum at 
the onset of the phase transition or the critical temperature following Borrill and Gleiser . 
Taking the initial condition as 

$(0,3;) =0, $(0,0;) = 0, (2.16) 

for all X , we follow the evolution of the field to trace the fraction of the symmetric phase, 
/o(t), which is defined by the fractional volume of the lattice space with $ < while 
the fraction of the asymmetric phase, /+(t) = 1 — /o(t), is that with $ > As in the 
non-selfinteracting case, the system is discretized and the second-order staggered leapfrog 
method is used. The noises are also generated in the same way. We have confirmed in all 
the cases of our interest the change of rj affects only the relaxation time scale but properties 
of the final configuration are insensitive to it |0. So hereafter we report the results with 
rj = 1. Contrary to Borrill and Gleiser §], we perform numerical calculations with various 
values of the grid spacing, which turn out to affect the final configuration greatly as seen 
below. 



First, as often assumed in the literatures 0] ||T^, we adopt so-called "correlation length" 
or the curvature scale of the potential at $ = 0, r = (f/"[$])~^/2, as the coarse-graining 
scale, namely, the lattice spacing. We thus set 5x = {9l - l)-^'^ and arrange a lattice with 
N = 64^, 6t = 0.1, and t = 1500 . For several A's, the fraction of the symmetric phase, /o(^) 
is depicted in Fig. ^. For A = 0.06 {6x ~ 7.9) corresponding to the Higgs mass Mh — 60GeV 
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TT|, the phase mixing does not occur. The result shows that if the Higgs mass is not too 



large, the phase mixing does not occur, which is consistent with the results of |14]. In 



Fig. ^ we have depicted the correlation function, ($(a;) $(?/)), at t = 1500, which does not 
necessarily approach zero for large separation due to the fact that the average value of $ 
is not equal to zero. In this figure we have also depicted theoretical curves Eq. ( p. 51) with 
m =([/"[$ = 0])^/^ in a dimensionless unit, namely. 



( ) = E / ,,,, ( m^{r^ + 2Dn^/9l) ) . (2.17) 



The former damps much more mildly than the analytic counterpart ( p. 17 ), motivating us to 
study the case with the smaller lattice spacing as has been done in [|]. 

Next we investigate the dependence of results on the grid spacing 5x. Two lattices are 
arranged, one with N = 64^, 6t = 1.0, t = 1500 and, 5x = 1.0, and the other with the 
same properties except for 6x = 0.5. The results are depicted in Figs. 5(a) and 5(b). The 
former case reproduces Borrill and Gleiser's result. Comparing both results, we find that 
the smaller the lattice spacing becomes, the phase mixing occurs for the smaller values of A 
corresponding to the lighter Higgs mass. This result can be understood as follows. Taking 
lattices is equivalent to cutting off the momentum. As is seen in Eq. ( p.4|) , the smaller 
lattice spacing we take, the larger momentum dominates and the more easily the phase 
mixing occurs. The correlation function are depicted in Figs. 6(a) and 6(b). In the former 
case with 6x = 1.0, it has a similar curve to the analytic one ( |2.17| ) except for the offset, 
while in the latter with 6x = 0.5 the correlation damps more rapidly. 

In order to examine dependence of the lattice spacing further, we have also performed 
numerical simulations with N = 64^, 6t = 0.1, t = 3000, and A fixed at 0.06. As is seen in Fig. 
1^, the results depend on the lattice spacings very much. 6x ^ 6.0 is critical and if we take 6x 
smaller, the phase mixing is manifest. Therefore unless we specify the lattice spacing from a 
physical argument, we cannot draw any quantitative conclusion about whether two phases 
are mixed or not. From the analogy with the non-selfinteracting massive scalar field analyzed 
in the beginning of this section, many people have adopted the Compton wavelength or the 
inverse-mass scale at $ = as the coarse-graining scale. In the present case, however, the 
correlation function changes significantly depending on the lattice spacing as shown in Figs. 
6(a) and 6(b) contrary to the case of non-selfinteracting field. Thus the above results suggest 
that the previous analyses adopting the correlation length, or the inverse-mass at $ = 0, 
as the coarse graining scale are inappropriate in the case of interacting potential. Thus we 
must reconsider the derivation of the Langevin equation before performing further numerical 



analysis based on the simplified equation (2.6) 



III. PROPERTIES OF THE NOISES DERIVED FROM A NON-EQUILIBRIUM 

QUANTUM FIELD THEORY 



Here we review a field theoretic approach to derive an effective Langevin equation with 
particular emphasis on the origin of its noise term. The standard quantum field theory. 
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which is appropriate for evaluating the transition amphtude from an 'in' state to an 'out' 
state for some field operator O, {out\ O \in), is not suitable to trace time evolution of an 
expectation value in a non-equilibrium system. In order to follow the time development of 
the expectation value of some fields, it is necessary to establish an appropriate extension of 
the quantum field theory, which is often called the in-in formalism. This was first done by 
Schwinger |jl6| and developed by Keldysh [T^. Here, following Morikawa |]l8l and Gleiser 



and Ramos [0, we first review briefly the derivation of an effective Langevin-like equation 
for a coarse-grained field using the non-equilibrium quantum field theory based on the in- 
in formalism and then extract necessary information on the noises which are essential for 
generating inhomogeneity of the system. 



A. Non-equilibrium quantum field theory 



Let us consider the following Lagrangian density of a singlet scalar field (p interacting 
with another scalar field x ^-nd a fermion ip for illustration. 



C 



1 2 2 



1 



1 



(3.1) 



Although the above Lagrangian density is much simpler than the standard model, it turns 
out that this model fully accounts the nature of bosonic noises arising from interactions 
with gauge particles and Higgs self-interactions as well as fermionic noises from quarks and 
leptons. 

In order to follow the time development of only the initial condition is fixed, and so 
the time contour in a generating functional starting from the infinite past must run to the 
infinite future without fixing the final condition and come back to the infinite past again. 
The generating functional is thus given by 



Z[J,V,V] = Tr 
= Tr 



T l^exp 
T_ ( exp 



P 



i I {J(p + Kx + r]il) + riil) 

i J {J-'-P- + K_x- + V-i'- + V-i'-] 



exp 



i / ( J+^+ + K+X+ + ^+^+ + 



P 



(3.2) 



where the suffix c represents the closed time contour of integration and v9+(x+, V'+j '^+) ^ 
field component ipi^x^i^^i^) on the +-branch (— cxo to +oo), (y9_(x-, V"-; V"-) that on the — 
branch (+oo to —oo). The symbol T represents the time ordering according to the closed 
time contour, the ordinary time ordering, and T_ the anti-time ordering. J, K, t], and 
fj imply the external fields for the scalar and the Dirac fields, respectively. In fact, each 
external field J+(i^+, ?7+, 7)+) and J-{K_,r]_,fi_.) is identical, but for technical reason we 
treat them different and set J+ = J_{K+ = K_,rj^ = ri_,f]+ = fj_) only at the end of 
calculation, p is the initial density matrix. Strictly speaking, we need couple the time 



7 



development of the expectation value of the field with that of the density matrix, which 
is practically impossible. Accordingly we assume that deviation from the equilibrium is 
small and use the density matrix corresponding to the finite temperature state. Then the 
generating functional is described by the path integral as 

Z[J,K,ri,f]] = exp(^iW[J,K,ri,f]]^ = J^V(fi J^Vx J^Vi; J^Vip* expiS[(fi,x,i^,^, J, K,ri,f]] 

(3.3) 

where the classical action S is given by 

S[(p, Xj'^j Jj K,ri,f]] = j d^x |£ + J{x)ip{x) + K{x)x{,x) + ri{x)^p{x) + f]{x)'ilj{x)^ . 

(3.4) 

As with the Euclidean-time formulation, the scalar field is still periodic and the Dirac field 
anti-periodic along the imaginary direction, now with (p{t,x) = ip{t — il3,x), x(t, a;) = 
x(t — i/3, x), and ^p{t, x) = —^p{t — i(3, x) [pO[] . 

The effective action for the scalar field is defined by the connected generating functional 

as 

r[0] = W[J,K,7],f]]- Jj^xJ{x)(f){x) , (3.5) 

where 0(a;) = 6W[J, K,ri,f]]/6J{x). 

We give the finite temperature propagator before the perturbative expansion. For the 
closed path, the scalar propagator has four components. 

'""'^yc-ix-x') G^ix - x') ) 

^ / Tr[T+x(x)x(a;')p] Tt[x{.x')x{,x) p] \ 
~ \ Tr[x(a;)x(a;')p] Tr[T_x(x)x(a;')p] / 

ik{x-x') ( G^{k) \ . . 



where 



^ k^ — mi + ie A\ / \ K' 

Giik) = 'I . + 2nn^ik) S{k' - m^) , 
^ — mi — le ^ 

G+{k) = 2Ti\e{-k,) + n^{k) ] 5{e - ml) , 

G-(k) = 2n[e{ko) + %(fc) ] - ml) , (3.7) 



with n^(fc) = (e^'^^(^) - l)-\ u^{k)^ = k^ + ml, and e{ko) = e{ko) - e{-ko) 0. Similar 
formulae apply for if field as well. 



8 



Also, for a Dirac fermion we find 



where 



/ S^{x — x') S^^{x — x') \ 
\ S^{x — x') S^{x — x') J 

I Tr[T^ip{x)'ijj{x')p] Tr[—'ijj{x')'ijj{x)p] 
y Ti['ijj{x)'ijj{x')p] ^T[T_^jJ{x)'ljJ{x')p] 

-ik(x-x') ( Si>{^) ^iiif^) \ 



(3.8) 



Slik) = — - 2irn^{k){/c. + m^) S{P - ml) , 

^ /t — + te ^ 

S^{k) = ^- : 2Txn^{k){/i + m^)6{k^ - ml) , 

^ /c — m^ — te ^ 

S^ik) = 27T[9{-ko) - n4k) ] + m^) - m\) , 

S^^ik) = 27r[e{ko) - n^k) ] + m^) 6{e - m^) , (3.9) 

with n^{k) = (e^^'^(^) + l)-\ uj^{k)^ = fc^ + |l . 



B. One-loop finite temperature eff"ective action 

The perturbative loop expansion for the effective action F can be obtained by transform- 
ing —>■ ifo + ( where ipo is the field configuration which extremizes the classical action 
S[(p, J] and C is small perturbation around ipo. Up to one loop order and 0{X^, g^, p), T is 
made up of the graphs as depicted in Fig. ||. Summing up these graphs, the effective action 
F becomes 



d^x |0A(a;)[-n - - M^](l),{x) - ^ (AMx)^lix) + 0c(x)0i(x) 



(j)A{x)(j)c{x)(j)l{x') + -(j)A{x)(j)c{x)(j)\{x') 



d'^x J d'^x' {Ai{x- x') + As{x- x') } 
-2 J d^x J d*x'A2{x-x')(f)A{x)(f)c{x') 
+- J d'^x J d'^x' {Bi{x — x') + B'j,{x — x') } 0A(a;)0A(a;')0c(a^)0c(a^') 

+B2{x - x')(f)A(x)(f)Aix' 



(3.10) 



where 



+ 



(3.11) 
(3.12) 
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(3.13) 




IvI 


X r l + 2n^(q) 
~ 2J (27r)3 2uo^{q) ' 




Ai{x 


-x') 


= ^lm[G^{x-x')^]e{t-t') . 
2 


(3.15) 


A2{X 


- X) 


= f^lm [ S^fjix -x)S^ [x - x)]0{t -t) . 


(3.16) 


Asix 


-x') 


\2 

= yIm[Gj(x-x')']^^(t-t')- 


(3.17) 


Bi{x 


-x') 


= |-Re[Gj(x-x'f]. 


(3.18) 


B2ix 


— x) 


= ~fRe[S^^{x-x')S^/{x'-x)]. 


(3.19) 


B3{X 


-x') 


= yRe[Gj(x-x')^], 


(3.20) 



The last term of ( p.lO| ) gives the imaginary contribution to the effective action T. We 
can attribute these imaginary terms to the functional integrals over Gaussian fluctuations 
and ^2 [|18l- That is to say, we can interpret that the imaginary part of the effective action 
comes from random fluctuations onto the expectation value. Thus we rewrite (|3.10| ) as 



exp(zr[0e,0A]) =/^^^l/ ^^6AN^2[6]exp{^S'eff[0c,0A,ei,6]} 



where 



S, 



cS\<Pc, < 



0A,^i,6] = Rer + / d^x[^i{x)(f)c{x)(f)A{x) + ^2ix)(f)Aix)] 



with the probability distribution functional 

1 



Mi exp 



d^x / d'^x'^i{x)B^ ^{x — x')^i{x') 



1,2). 



(3.21) 
(3.22) 

(3.23) 



where A/^ is a normalization factor and Bi{x — x') = Bi{x — x') + B^lx — x'), B2{x — x') 
B2(x-x'). 



C. Equation of motion 

Applying the variational principle to Ses, we obtain the equation of motion for 



6 4) A 



0. 



<Pa=0 



(3.24) 



^Prom ( g:22D and (|3l0| ), it reads 



A 



a + M^ + M^)<j),{x) + -(l)lix) 



+ (j)c{x) I d^x' I dt'Ai{x - x')(f)l{x') + 2 / d^x' / dt'A2{x - x')(f)^{x') 

= 0e(x)ei(x)+6(x), (3.25) 
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and 



{Ux)U^')) = - x') , (3.26) 

where Ai{x — x') = Ai{x — x') + A^lx — x') . Though Ai and Bi has two contributions from x 
and if fields, they have the same properties except for the values of coefficients and masses. 
From now on, we consider only the contribution from x field for simplicity and omit the 
suffix c. The right-hand-side of ( |3.25|) are the noise terms, while the last two terms of the 
left-hand-side are combination of a dissipation term and one-loop correction to the classical 
equation of motion which would reduce to a derivative of the effective potential, V^g^cp), if 
we restricted to be a constant in space and time. 

The above equation (|3.25| ) is an extension of equation (3.2) of Gleiser and Ramos []I9 



m 



that we have incorporated not only self-interaction but also interactions with a boson x and 
a fermion ip. In Gleiser and Ramos proposed to adopt several further approximations 
to reduce their equation to the form of the simple Langevin equation like ( |1 . 1|) and ( [1.21) . 



In particular, for the purpose of simplifying the equation to a local form they handled 
spatial nonlocality by considering only contributions with zero external momentum, which 
is physically equivalent to dealing only with nearly spatially homogeneous fields. With this 
approximation the correlation function of the bosonic noise (|3.18|), for example, becomes 



(g;e*<"-"'' / ||Re |GJ(,, t - t')G^^(g -k.t- t') 



,4 



k=0 ' 

^6%x-x')I^^Re[G^{q,t-t')r. (3.27) 

We thus obtain spatially uncorrelated noise, which would violate spatial homogeneity of 
in the severest manner and lead to a self-inconsistent result. 

Since the noise term in ( |3.25|) are the only source of inhomogeneous evolution of (p, we 



should not adopt the above approximation ( |3.27| ) for their correlations. Instead, keeping the 
original form of the correlation functions of the noises, ( p.26| ) with ( |3.18| ) through ( p. 20]) , 



we can obtain important informations about inhomogeneous fiuctuations which help us to 
set the lattice spacing of simulations. For this purpose we calculate the spatial correlation 
functions of the noises more explicitly in the next subsection. 



D. Spatial correlation of noises 

/^From the above discussion, we see the correlation length of the noise is the most im- 
portant scale in order to investigate the effect of fiuctuations on the dynamics of the phase 
transition. The correlation function of the noises are given in ( p.26| ), which are, unfortu- 
nately, too complicated to apply for numerical simulations directly. Since nontrivial tempo- 
ral correlation is expected to affect the relaxation process and its time scale only, we adopt 
an approximation that temporal correlation is white but take spatial correlation fully into 
account. So we evaluate the equal-time spatial correlation. 

First we consider the case of the bosonic noise. The equal-time propagator in momentum 
space is given by ^| 
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Gj(fc,0) 



1 



2uJk 
1 

2uju 



1 + 2E 



n=l 



Then that in configuration space propagator reads 



(27r) 

(27r)3 2u;fc 



1 + 2E 



n=l 

oo 



iriKr)+2E 



(3.28) 



, (3.29) 



which of course has the same form as (|2.5| ). Using the above representation, the equal-time 
spatial correlation of the bosonic noise becomes 



;ei(a^)6(0))e,.al- 



9^ 



time 



Re 



2 4 r 



Ki {m^ Jr"^ + 11? 13"^ 



(3.30) 



which damps exponentially for r > m^^ . 

A, 

For the massless bosonic noise, we can calculate the sum of the infinite series to yield 



Gj(a;,0) 



coth ( -vr 1 . 



(3.31) 



and the equal-time spatial correlation of the bosonic noise becomes 



:ei(^)ei(o) 



■4 1 



coth -vr 



9 



327r2/52 r2 



for , r 3> — . 

TT 



(3.32) 
(3.33) 



We thus find that for the massive bosonic noise the correlation function damps exponentially, 
particularly the damping scale is the inverse of mass, while for the massless bosonic noise it 
damps much less rapidly, according to a power-law. 

For the fermionic noise, similarly, the equal-time propagator in momentum space is given 

by m 



<(fc,o) 



1 



2ujk 
1 

2uJk 



(m^ - 7 ■ fc) 
(m^ - 7 • fc) 



i + 2E(-i 



n=l 



(3.34) 
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and that in configuration space reads 



,0)=/ 



l + 2^(-ire 



n -nl3wk 



n=l 



im: 



47r2 



r r 



ii'2(mv,r) + 2^(-ir 



n=l 



(3.35) 



Also, the equal-time correlation of the fermionic noise is given by 



;6(a^)6(0))equal- 



time 



-f Re 



Trfe(a^,0)5^(-a;,0) 



i^2(mv,r)+2^(-ir 



n=l 



47rV2 



n 2 



n=l 



(3.36) 



which damps exponentially at the inverse mass scale for m^P ^ 1, and at the scale (3 for 
m^(3 < 1. For the massless fermion, we find 



7 ■ cc 



TT 



+ 



TT" cosh (^^tt) 



47r2r2 [ /3r sinh (|7r) sinh^ (^tt) 



(3.37) 



and 



■time 



47^4^2 



TT 



TT" COSH ( -TT 



/?'^sinh(|7r) /?^sinh2(^7r 



f e-- 

4^4 ^2 



tor , r ^ — . 

TT 



(3.38) 



(3.39) 



Note that for the fermionic noise, unlike the bosonic case, the correlation function damps 
exponentially regardless of the mass. This is the very interesting feature of the noise. We can 
physically interpret this feature as follows. Since there is Pauli's blocking law for a fermion, 
particles arc apt to separate from one another and the correlation is easily destructed. On the 
other hand, bosonic particles can occupy the same state and the correlation is kept. When 
the temperature is zero, both fermionic and bosonic massless propagators damp according to 
a power-law. The above feature is an example of the fact that statistical difference between 
fermions and bosons appear more markedly at finite temperature than at zero one. 
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E. Dissipation term 



The equation of motion ( |3.25| ) derived in the subsection |III C| has contributions repre- 
senting the dissipative effect in the last two terms of the left hand side. Since these terms 
are nonlocal in time, what is often done in the literatures [0 [1^ to extract local terms 
proportional to is to assume that the field changes adiabatically, or put 



(f{x', t') ^ (P'\x\ t) - n{t' - t)(P''-\x\ t)^{x\ t) 



(3.40) 



in the integrand of (|3.25| ). But the dissipation terms thus evaluated vanish as long as 
we use the bare propagators. This is usually interpreted as a manifestation of the fact 
that the dissipative effect is intrinsically a non-perturbative phenomenon and cannot be 
investigated from the perturbation theory [22|. In order to see a damping effect, we must 



observe the system for a finite duration of time typically proportional to the inverse of some 
coupling constants. After this period, however, the perturbation theory may have broken 



down as explained in |^ using a toy model. So, in order to find the damping effect, non- 
perturbative terms are often incorporated by using a "dressed" propagator with an explicit 
width, which is obtained by resuming higher-loop graphs of some classes, instead of the bare 
one |[T8| , p^ . Although we may obtain finite dissipation terms in this way, they do not yet 
satisfy the fluctuation-dissipation relation. More serious is the fact that these approach is 
not self-consistent as criticized by Gleiner and Miiller p3|. In fact, apart from the validity 
of perturbation, the adiabatic expansion ( |3.4CI| ) itself breaks down before we can observe 
dissipation effect |^ 



In order to cure the situation, Gleiner and Miiller proposed to adopt the linear 
harmonic expansion of the Fourier mode (j){k,t') as 



{k,t') ~0(fc,t) cos[iUkit' -t) 



[k,t) — sin [ij kit' -t)] 



(3.41) 



and calculated the dissipation term in a simple model. They have shown that the fluctuation- 
dissipation relation is met in the classical limit |^ . Unfortunately, however, the applicability 



of this method is rather limited and we cannot calculate the dissipation term in our model 
explicitly. Therefore, here we make much of the thermodynamics and derive it so that the 
fluctuation-dissipation relation is met. We also identify the remaining part of the integral 
terms of ( p.25| ) with the derivative of the effective potential V^fj(0) , primarily for simplicity. 
But if the system turns out to be homogeneous as a result of numerical simulations, this 
choice will be justifled since the homogeneous expectation value should take a value where 
the effective potential is minimized, but otherwise unjustifled. In the latter case the one- loop 
approximation also breaks down and then we would have to deal with the full effective action, 
which is beyond the scope of the present analysis. We thus interpret the integral terms in 
( j.25 ) as consisting of two parts, one contribution to the derivative of the effective potential 
and the other to the dissipation term, and this division is done so that the dissipation 
term satisfles the fluctuation-dissipation relation. Note that, although the above procedure 
is physically motivated, it should be regarded as an ansatz rather than an approximation 
and we have been unable to derive these terms rigorously from flrst principles at present. 



14 



Nevertheless we stress that our approach is self-consistent if the field configuration turns out 
to be homogeneous. 

We next put the above scheme into practice. First we consider the case thermalization 
proceeds only through bosonic noises and determine the corresponding dissipation term. To 
do this we rewrite ( p.25| ) in the form 

□(/)(a;) + Kfr(0) + <l>ix) J d'^x'Aiix - x' , t)(j){x', t)<j){x' , t) = (j){x)^i{x) , (3.42) 

where Ai{x — x', t) is to be determined so that the fluctuation-dissipation relation is met. For 
the purpose of applying to numerical simulations we adopt an approximation that ^i(x) is a 
temporally white noise with ( ^i(a;)^i(x) ) = Bi{x—x')5{t—t'), which is a good approximation 
since temporal correlation damps exponentially beyond \t—t'\ > /?/(27r). Then this Langevin 
equation can be converted into the Fokker-Plank equation. 



dW 



d-^x 



6n 6W 5H m 

+ 



+ 



b-Kix) 6(f){x) 6(f){x) 6tt{x) 
6 



57r(x) 



x,t) J d^x'Ai{x - x',t)(f){x',t)n{x',t) -W 



+ 



^Jd^x J d^x'4){x,t)4){x',t)Bi{x - x'] 



6^W 



5tt{x)6tt{x') 



(3.43) 



where W[(j){x, t), n{x, t)] is the distribution function, tt{x) = 4>{x), and H is the Hamiltonian, 
which is given hj H = J d^x | vr^ + |(V0)^ -|- Kfr(0) • In order that this equation has a 
stationary solution, it is at least necessary that Ai{x — x',t) does not depend on time, 
namely, Ai{x — x', t) = Ai{x — x') . We require that W^t = Mt^~ constitutes a stationary 
solution of the above equation. 



= \ d^x \ d'^x'(j){x,t)(f){x',t) 



6^{x - x')W - -n{x)TT{x')W 



Ai{x - x') - 



Then we find 



Ai(x - x') 



1 

2T 



Biix-x') 



and 



2T 



a^{x) + VM) + t;;^^^) / d^x'Bi{x - x')<p{x' ,t)^{x\t) = ^{x)^i{x) 



The fermionic contribution can be treated similarly, and we find 
□0(x) + Kff(0) + ^/ d'x'B^ix - x')^{x',t) = 6(x) 



Bi{x-x')' 
2T 

(3.44) 
(3.45) 

(3.46) 
(3.47) 



{U^)U^)) = B2{x-x')5{t-t'), 
in the case thermalization is realized only through fermionic noises. 



(3.48) 
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IV. NUMERICAL SIMULATIONS 



Using the phenomenological Langevin equation obtained in the previous section, we now 
perform numerical simulations in the same way as in Sec. |I[ The essential difference from 
Borrill and Gleiser's formulation is that we generate noises which have spatial correlations 
calculated in the previous section. 

In order to approach some aspects of the electroweak phase transition in the standard 
model, we adopt the one-loop improved effective potential of the Higgs field in our phe- 
nomenological Langevin equation. Since the properties of massless bosonic noises such as 
those from gauge interactions and fermionic noises are very different from each other, we 
analyze the cases thermalization proceed through bosonic noises and fermionic noises sepa- 
rately. 



A. Bosonic noise 

We first examine an extreme case that thermalization of (j) proceeds by virtue of only 
massless bosonic noises which has a power-law correlation function such as those arising 
from gauge interactions. As a power-law is scale free, we set the same lattice spacing as the 
fermionic case, /?/(27r), in order to see the difference between the behavior of bosonic noises 
and that of fermionic counterparts. 

Before calculation it should be noted that the bosonic noise we have derived so far is 
multiplicative and accordingly, if we set the initial condition as $ = $ = 0, the system does 
not evolve. Hence, in this case we add the following contributions with two loops and order 
g'^ (Fig. which lead the additive noise. Then the correction to the effective potential is 
given by 

Ar[0c, 0a] = — / d'^x(f)/\{x)/^V4)c{x) — 2 / d'^x / d'^x' Aj^{x — x)(j)/\{x)(j)c{x) 



where 



I 



-\-- I d^x I d'^x' Bi{x — x')(j)/\{x)(j)/^{x ) , (4.1) 



^V='4[dt'[ T^Im \ Glik, t-t')f Bit - t') [ > (4.2) 



2 J J (27r)3 L XV ' n ^ 'J (27r)3 2uo^{q 



A^{x - x') = y Im [ Gj(x - x'f G^ix - x') ] e{t - t') , (4.3) 

B^{x-x') = ^Re[Gl{x - x'f Gl{x - x')] . (4.4) 
2 

The last term is imaginary and we regard it as the term coming from the stochastic noise as 
before. The dissipation term is obtained so that the fluctuation-dissipation relation is met. 
Then the dimensionless equation of motion only with bosonic contributions becomes. 
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[X) 



V2$(x) - 



dU{^) 



~^x) [ d^x'Biix - x')^x',t)^x',t) - 4 / d^x'Bi{x - x')^x',t) 

+^x)U^)+U^), (4.5) 



where 



) = B,{x - x')S{t - t') = T-—jr-,S{t - t') , 



X — X' 



iU^M^')) =B^{x-x')6{t-t') 



Or 



AttIx — X 



(4.6) 
(4.7) 



In reality, the correlation function of the noises as derived from a fundamental theory is 
proportional to some powers of coupling constants. However, once the fluctuation-dissipation 
relation is assumed, its magnitude does not affect the final equilibrium configuration at all. 
Hence we have normalized the amplitude of noises as above in our numerical calculations. 

Noises with the above correlation can easily be generated working in Fourier space. Since 
we are assuming ^{x) is random Gaussian, its Fourier transform. 



ak)^jd'xe-^^-''i{x), 



(4.8) 



also satisfies the Gaussian probability function, and the distribution function for each mode 
becomes 



P[C(fe)] = Ar'exp[-i| 



d'k |C(fc)| 



2 n 



(27r)3 P{k) 

Here A^' is the normalization factor, and P{k) is the power spectrum given by 

P{k)^ j d^xe-'^-^B{x) , 



(4.9) 



(4.10) 



and we have used the reality condition for ^{x). This distribution has no correlation between 
each Fourier mode. Therefore we have only to generate the white noise in the momentum 
space and Fourier transform it into the configuration space. 

Next we discretize the system. Then the discretized master equation becomes 



l,n+l/2 



+ 5t V'$,- „ - 



dU 



Sx 



3 
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s=x,y,z {Sx) 

where we have used the second-order leapfrog method and adopted the Crank-Nicholson 
scheme only for the diagonal term because of its dominance for numerical stability. We also 
set B^^ = = except B^q = 1/N'^,B^q = Initial conditions are given 

by $^ g = $^ Q = for all i's. 

The results are depicted in Fig. Massless bosonic noises do not disturb the homoge- 
neous field configuration at least for small enough values of A. 



B. Fermionic noise 



Next we consider the fermionic case. The spatial correlation of the fermionic noise 
damps exponentially for both massive and massless fermions. The fermions interacting with 
the Higgs field in the standard model acquire finite masses if and only if the Higgs field 
has a non-vanishing expectation value. Hence in the same spirit as adopting the effective 
potential in the equation of motion of the field, we perform numerical calculations with 
massless fermionic noises. That is, if remains to be equal to zero homogeneously reflecting 
the initial condition, this choice is consistent. In fact, however, even when phase mixing 
is manifest, the expectation value of the Higgs field remains at most about 50GeV for 
Mh — 60GeV, which means that quark masses remain much smaller than the temperature 
and massless approximation itself is justified in this well. 

Since the correlation of the fermionic noise damps exponentially, we can use the white 
noise as long as we take the lattice spacing larger than the correlation length, which is equal 
to /9/(27r) and in our dimensionless unit it corresponds to Sx = (2D)^/^/(27r6'c) (= 0.092 for 
Mh = 60GeV). Then the master equation is the same as that used by Borrill and Gleiser 
or equation ( p.l4|) . 



where 

iU^M^') ) = "^vO^it - t')6^{x - x') . (4.13) 
The results of the solutions of this master equation has already been depicted in Fig. |^ 



for various grid spacing 6x. Figure |TT| represents the case 6x is set to be the fundamental 
length. As is seen in these figures, the fermionic noises are more effective to disturb the field 
configuration from a homogeneous state to an inhomogeneous one with possible mixing of 
two phases. 
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V. ANALYTIC INTERPRETATION OF NUMERICAL RESULTS 



As in Fig. 0, in some choice of 5x and A the system seems to relax to a state with /o 
some value between /q = 0.5 and 1. One may wonder that this is because our simulation 
time is so short that the system still keeps the memory of the particular initial condition, 
and that if we could trace the evolution for a long enough period the system would relax 
into a state with /o = 0.5. In fact, however, if we examine the time variation of the field 
configuration by observing its snapshot at different times as depicted in Figs. 12(a) - (e), we 
can easily convince ourselves that the system is in a stationary state with a constant /o(> 0.5) 
repeating creation and annihilation of a number of small domains of the asymmetric state 
whose typical radius is at most a few times the lattice spacing. In this section we would 
like to present an analytic argument to support that the state we followed in the numerical 
simulation is a thermal state. Similar analysis has already been done by Gleiser, Heckler, 



and Kolb in a slightly different situation. See also Gelmini and Gleiser |]25[ . 

Let (i+(-R, $,t) be the number density of the asymmetric-state bubbles with radius 
R{> 5x) and amplitude $(> $_) at t. In order to obtain the Boltzmann equation for 
$,t), we count the processes which change it: i) Thermal nucleation of asymmetric- 
state bubbles in an almost homogeneous symmetric-state sea. ii) Annihilation of these 
asymmetric-state bubbles into the symmetric-state sea. This is to be distinguished from 
the process of nucleation of a symmetric-state bubble in a homogeneous background of the 
asymmetric-state whose rate would be identical to that of i) in a degenerate potential. We 
expect the process i) has smaller rate than the process ii), since the former requires more 
energy. We do not take into account the process that nucleated bubbles dynamically shrink 
{\v\dd^/ dR term in p4|) because the typical radius of nucleated bubbles is comparable to 
the lattice spacing and shrinking bubbles and vanishing bubbles are hardly distinguishable 
in our simulations. Then the Boltzmann equation for d^{R, $,t) becomes 

dd^{R,^,t) 4:7iR^ 

= (1 ~ f+it)) G^(s-phase^R,$) ' d+{R, $, t) G'(s-phase<S=R,*) , (5.1) 

where G(s-phase=>R,<i') is the nucleation rate for the process i), G'(s_phase<;=R,<i') is that for the 
process ii). We also assume that nucleation rates can be obtained from the Gibbs dis- 
tribution, namely G = Aexp{—F/9c), where A is a constant. For G (s- phase ^r,$) we put 
F = b^^R taking surface tension of the created bubbles into account with b a constant. 
Since the inverse process is not Boltzmann suppressed, we assume its rate is a constant: 

G(s-phase<^R,4>) = B. 

In an equilibrium state, dd+/dt equals zero for all $'s and i?'s, 

(1 - /f ) ^(-P*--^^'^) = ^rf_^(i?,$,t) . (5.2) 

^(s— phase-4=R,<I>) "J 



Summing these equations for all $'s and -R's leads to 

d^dR =11 — 

ISx J<S>- ' ' ' (j^(s-phase<^R,$) 



^(s-Phase^R.^) ^^^^^ r r^d^^R^<^^t^dMR = fir (5.3) 
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Then becomes 



where 



oo /-oo 



5x J<3?_ G(; 

/•oo 

/ exp{-F/9,)dMR. (5.5) 

JcJ>„ 



(s— phase<^R,'I>) 

oo roo 



is depicted in Fig. |1^ as a function of A. Since the case for 6x = 1.0 is suitable for seeing 
the change of the asymmetric fraction from non- mixing to percolation, we set 5x = 1.0. Then 
we can fit the analytic solution with the numerical simulation very well with A/B = 65 and 
b = 2.77. 

Note that even if |f | equals zero, f!^ can become values different from 0.5. The essence 
lies in the fact that creation rate of an asymmetric-state domain in the symmetric phase and 
that of its reverse can be different due to the surface tension even at the critical temperature 
if the background is sufficiently homogeneous. Of course, in the case the average amplitude 
of fluctuations is large enough, percolation occurs quickly and the two processes will have 
the same rate, resulting in = 0.5. 

The above argument is expected to apply only when $ is localized around the origin 
initially. If it is localized around $ = initially, on the contrary, we expect that the 
system will settle into a state with = 1 — //(/ + 1) or /q"^ = //(/ + 1), since the potential 
we are using is symmetric. In order to see what happens in the case $ is not localized 
around either minima, we have run five simulations starting from a checkerboard made of 
$ = and $ = using different realization of random numbers. The result is depicted in 
Fig. [1^ for A = 0.06, 6x = 7.0 and = 64^, which shows that the system approaches to 
either equilibrium state with or /g'^ = //(/ + 1), although it takes much longer time to 
relax than in the cases with $ = or $ = initially. 

This result implies that the configuration with /+ = /q = 0.5 is unstable in this case 
even if it contains maximum number of microscopic states. Note that this configuration 
also costs more energy of domain boundaries than any other configuration. This is why the 
system may relax to a configuration with /+ 7^ /q. 

Although these arguments are interesting in themselves, we must be cautious with their 
interpretation, that is, it may not be directly relevant to the actual dynamics of EWPT 
because our use of the effective potential in the Langevin equation is not strictly justifiable 
in the case field configuration becomes inhomogeneous. The same warning also applies 



to Anderson's argument p8[, who claims that basic picture of phase transition through 
subcritical bubbles is in contradiction with the second law of thermodynamics. This criticism 
is also based on an expression of the free energy which is not strictly correct in inhomogeneous 
situations. 



VI. SUMMARY AND DISCUSSION 
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In the present paper we have performed a series of numerical simulations of the Langevin 
equations toward understanding aspects of a weakly first-order cosmological phase transition 
such as the electroweak phase transition. 

First we have confirmed that the simple Langevin equation (|1 . 1| ) with random Gaussian 
white noise ( |1.2| ) can reproduce thermal equilibrium state of a massive non-selfinteracting 
scalar field, in particular, that the correlation length is given by the inverse-mass scale 
independent of the lattice spacing. We have then applied the same technique to one-loop 
improved effective potential of the Higgs field in the electroweak theory. Taking the coarse- 
graining scale or the lattice spacing equal to the inverse-mass scale at = 0, we have 
confirmed that phase mixing does not occur at the critical temperature for a small enough 
Higgs mass, consistent with the previous analytic estimate of the amplitude of fluctuations 
[|I^ [|T5|. At the same time, we have argued that in order to reproduce the shape of the 
corresponding massive scalar correlation function we should take 5x smaller. As a result 
of such simulation we have found that the correlation function of the Higgs field obtained 
numerically may damp at smaller scale depending on the choice of 5x, so it has proved that 
the so called "correlation length" or the inverse-mass scale at = is not a good measure of a 
coarse-graining scale in the case the potential contains a nontrivial interactions. In this case 
we have also found the final configuration of the simulation depends on the lattice spacing 
severely, which is a manifestation of the fact that lattice spacing serves as an ultraviolet 



cutoff of otherwise divergent theory. Since the renormalization prescription of ||2^ does not 
work to the temperature-dependent potential, we have tried to fix the lattice spacing from 
a physical argument. 

For this purpose we have reexamined derivation of the Langevin-like equation in the 
literatures |18| We stressed the importance of the correlation of the noise terms which 



are the only source of inhomogeneity in the system and so the simulation should be done 
reflecting their properties. Although the noise terms can be derived from the perturbative 
non-equilibrium field theory, no completely satisfactory derivation has been given of the 
other ingredient of thermalization, namely, the dissipation terms. Hence we made much of 
thermodynamics and determined them from the fluctuation-dissipation relation. Another 
difficulty of the equation of motion from the perturbative effective action is that it contains 
integral terms which are non-local in both space and time. Since it is impossible to deal 
with them numerically, we have replaced them with the derivative of the one-loop effective 
potential. This procedure may be justified if and only if the field configuration remains 
homogeneous. Since we set a homogeneous and static initial condition, 0(x, t) = 0(x, t) = 0, 
and what we are concerned is if the system remains homogeneous, the above approximation 
is sensible. Note, however, that in the case phase mixing is manifest in the final result, our 
simplified equation is no longer valid. Then the one-loop approximation also breaks down at 
the same time and we would have to deal with the full effective action, which is formidable 
now. 

Keeping the above-mentioned limits of our approach in mind, let us consider implication 
of the results of numerical calculations. We examined the effects of bosonic noises and 
fermionic noises separately. In the former case the field remained practically homogeneous 
at least for small enough values of Mh, and in the latter case phase mixing was evident 
and our approximation broke down. If only one species of noises and the corresponding 
dissipation term are taken into account, the strength of the noise and the dissipation do 
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not affect the final equilibrium configuration, as far as the fluctuation dissipation relation is 
satisfied, although they do affect the relaxation time scale. In the realistic case, both types of 
noises are present and their amplitude is also important to determine thermalization process 
of the system. Since the Yukawa coupling of the top quark is larger than the square gauge 
coupling and fermionic noises are more effective, we expect that the results of subsection [IV B| 
apply to the actual electroweak phase transition. In short, although gauge interaction plays 
the essential role to induce the cubic term in the effective potential, the non-equilibrium 
dynamics is dominated by fermionic interactions, and as a result, the conventional picture 
of first-order phase transition based on the one-loop potential is suspect. 

So far we have traced the behavior of the expectation value of the scalar field. Although 
the equation of motion has been motivated from quantum theory, as far as we concentrate 
on an expectation value, we must make sure that the effect of quantum uncertainty is 
sufficiently small. This condition is satisfied if the number of quanta contained in one lattice 
volume is much larger than unity The simulations with a lattice spacing equal 

to the fundamental correlation length of the fermionic noise, (3/{2tt), does not satisfy this 
constraint. Our conclusion, however, remain intact since at the critical value of 6x = 6.0 in 
Fig. ^, one lattice volume already contains more than 40 quanta, much larger than unity. 

In the light of the above results and from the fact that all the previous literatures claiming 
that usual nucleation picture in the homogeneous background applies to the electroweak 
phase transition have estimated fluctuation on too large a spatial scale, it is now evident 
that the conventional picture with one-loop effective potential does not work. On the other 
hand, in order to clarify real dynamics of the phase transition we must say that much work 
is to be done including derivation of correct equation of motion of the order parameter which 
should replace the crude phenomenological equation employed here. 
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FIG. 1. Correlation function {^{x)^{y)) for the massive non-selfinteracting field at finite tem- 
perature. 
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FIG. 3. Fraction of the symmetric phase, /o, with the lattice spacing 5x taken as the curvature 
scale of the potential at $ = 0. {5x = 7.2(A = 0.050), 7.9(A = 0.060), 8.6(A = 0.070)) 
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FIG. 4. Correlation function ($(£c)$(y)) at i = 1500 for the interacting scalar field with the 
lattice spacing 5x taken as the curvature scale of the potential at $ = 0. 
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FIG. 5(a) Fraction of the symmetric phase, /o, for 6x = 1.0. 
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FIG. 5(b) Fraction of the symmetric phase, /o, for 6x = 0.5. 
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FIG. 6(a) Correlation function ($(£c)$(y)) for the scalar field at with the one-loop electroweak 
potential {6x = 1.0). 
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FIG. 6(b) Correlation function ($(£c)$(y)) for the scalar field at with the one-loop electroweak 
potential {Sx = 0.5). 
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from the bottom to the top with A = 0.060. 
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FIG. 9. Contribution to F with two loops order and 0{X'^,g'^, /^). 
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FIG. 12(a) Snap shot in a plain at t = 1000 of the simulation with Sx = 1.0 and A = 0.02. 
Black dots represent asymmetric-state bubbles, while white ones a symmetric-state background. 
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FIG. 12(b) Snap shot at t = 2000. 
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FIG. 12(c) Snap shot at t = 3000. 
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FIG. 12(d) Snap shot at t = 4000. 
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FIG. 12(e) Snap shot at t = 5000. 



40 



.5 




Lambda 



FIG. 13. Comparison of analytical and numerical results. 
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FIG. 14. Fraction of the symmetric, /o starting from a checkerboard. 
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